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ABSTRACT 



o 

Context. Paper I presents the magnetic structure, inferred for the photosphere and the chromosphere, of a filament that developed in 
Cj-l! active region (AR) NOAA 10781, observed on 2005 July 3 and July 5. 

Aims. In this paper we complement those results with the velocities retrieved from Doppler shifts measured at the chromosphere and 
the photosphere in the AR filament area. 

Methods. The velocities and magnetic field parameters were inferred from full Stokes inversions of the photospheric Sii 10827 A 
line and the chromospheric Hei 10830 A triplet. Various inversion methods with different numbers of atmospheric components and 
different weighting schemes of the Stokes profiles were used. The velocities were calibrated on an absolute scale. 
Results. A ubiquitous chromospheric downfiow is found in the faculae surrounding the filament, with an average velocity of 
1.6kms _1 . The filament region, however, displays upflows in the photosphere on both days, when the linear polarization (which 
samples the transverse component of the fields) is given more weight in the inversions. The upfiow speeds of the transverse fields 
in the filament region average -0.15 km s~'. In the chromosphere, the situation is different for the two days of observation. On July 
O || 3, the chromospheric portion of the filament is moving upwards as a whole with a mean speed of -0.24 km s~'. However, on July 5 

only the section above an orphan penumbra shows localized upfiow patches, while the rest of the filament is dominated by the same 
downfiows observed elsewhere in the facular region. Photospheric supersonic downflows that last for tens of minutes are detected 
below the filament, close to the PIL. 

Conclusions. The observed velocity pattern in this AR filament strongly suggests a scenario where the transverse fields are mostly 
dominated by upflows. The filament flux rope is seen to be emerging at all places and both heights, with a few exceptions in the 
chromosphere. This happens within a surrounding facular region that displays a generalized downfiow in the chromosphere and lo- 
calized downflows of supersonic character at the photosphere. No large scale downfiow of transverse field lines is observed at the 
^ ' photosphere. 

Key words. Sun: filaments, prominences - Sun: faculae, plages - Sun: photosphere - Sun: chromosphere - Techniques: polarimetric 
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• 1 . Introduction shearing motions of the footpoints to form the flux rope in the 

corona with the aid of reconnection processes identified in the 
^ ! Velocity measurements in active region (AR) filaments are ex- photosphere as cancellation events. On the other hand, in the 
I ■ tremely scarce in the literature. Most studies have been carried FRE model the flux rope is formed deep in the convection zone 
• • | out only in quiescent (QS) filaments (outside ARs). However, (CZ) before it emerges through the surface and makes its way up 
_ ^ , quantifying velocities and plasma flows in filaments is a fun- into the corona. However, simulations based on buoyancy insta- 
damental step towards understanding their formation, evolution bilities encounter difficulties when lifting up a ma ss-loaded flux 
and disappearance. Filaments are structures formed from cool rope from the photosphere into t h e corona (e.g., iMurrav et"ail 
C3 ■ dense plasma that typically lies bet ween opposite polarities, i.e ., l2Q06h iMac Taggart & Hoodll20irj) . iKuckein et alj d2012l) (here^ 
at the polarity inversion line (PIL; Babcock & Bab cock|[l955l) . after Paper I) provide a more thorough discussion and additional 
in the chromosphere or corona. Due to the temperature contrast references to these t wo models (see also the extensive review by 
with their surroundings, they appear as dark structures when ob- iMackav et al.ll2010T) . 

served on the solar disk. Filaments are called prominences if ob- NT , , . , , .. . , , 

. . , r , Nowadays, high-resolution imaging has shown that qui- 

served in emission above the limb, but both terms are often used . c1 «. . c 11 1 r 

escent filaments are made of smaller-scale structures (e.g., 

interchangeably in the literature. bemouUn et alJ[l987t iLin et all I20051 and references therein). 

To date, there are two proposed scenarios for the forma- Doppler measurements in these quiesc ent filaments (e.g., 

tion process of filaments: the sheared arcade (SA) model and Mart res et al.lll98lt ISchmieder et a l. 2010) have revealed that a 

the flux rope emergence (FRE) model. Both result in very sim- wide range of velocities coexist in the sam e structure. For exam- 

ilar global structures, i.e., a flux rope that hangs above the pho- pie, the recent study bv lChae et alJ (120071) reported line-of-sight 

tosphere (although the SA model combines dipped field lines (LOS) velocities between +5 km s and +15 km s~' in the spine 

and flux ropes). However, it is in the formation process of the and in thr eads (respecti vely) of a QS filament (at coordinates 

flux ropes where the two models differ. The SA model requires N15-E26). lMeinl ( 1 19771) found that the axis of an AR filament 
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(at N7-W30) was at rest, while downflows on one side and up- 
flows on the other side were detected. Upwar d motions have also 
been d etected below QS and AR filaments bv lloshpa & Obridkol 
(1999). However, the main limitation of Doppler line shifts, and 
hence of the retrieved LOS velocities, is that the inferred val- 
ues are position-dependent all over the solar disk. Therefore, it 
is hard to compare Doppler velocities measured in different fila- 
ments. In addition, a thorough calibration, i.e., setting a correct 
zero value for the rest wavelength, is necessary, especially when 
calculating photospheric velocities (be low the filam ent) which 



are much smaller (< 0.3 km s 1 ; iMartres et al.l 



1976b than those 



Welsch & Fisher 



found in the chromosphere and corona (see 
120121) . 

It is presently unclear whether there is an unequivocal 
way of distinguishing which of the different filament for- 
mation scenarios applies to the observed AR filament. In 
this paper we propose to shed some light on this topic by 
studying the observed Doppler velocities inside the AR fila- 
ment that was analyzed in Paper I. In particular, we can ask 
what the expected velocity patterns (in and below the fila- 
ments) associated wi th these two models are. The flux rope 
models proposed by Iro n Ballegooijen & Martens ( 1989) and 
Ivan Ballegooijen & Mackavl d2007l) are built through successive 
photospheric flux cancellations driven by motions that force the 
field lines to converge towards the PIL. After the oppositely di- 
rected field lines are brought together, reconnection processes 
above the photosphere create a flux rope in the corona. The au- 
thors suggest that transverse magnetic field lines (loops which 
are almost perpendicular to the PIL) submerge at the PIL, and 
are transported downwards into the CZ as a result of the recon- 
ne ction processes mentioned before ( see panels d and e in Fig. 1 
of Ivan Ballegooijen & Martens|[T989l) . The typical size of these 
submerging loops is small (~ 900 km), which is large enough to 
be detected with currently available instruments. 

The prominence mod el presented in the work of 
iDeVore & Antiochosl (§000) does not show convergence of 
flux at the PIL. The helicity in the field lines is built up in the 
corona as a consequence of reconnection of the highly sheared 
field lines (twisted and stretched by the sheared displacements 
of the footpoints) with the overlying coronal arcade. The 
authors point out that the resulting magnetic structure is stable; 
hence, no upflows or downflows are expected in this model. 
However, the observed evolution of filaments suggests that they 
are of a non-steady nature. This calls for models with built-in 
evolutionary processes that allow the filament to disappear 
from the solar surface, transporting magnetic fields upwards or 
downwards through the photospheric boundary. 

The first complete description of the emergence of a flu x 
rope AR filament was presented by lOkamoto et al.l d2008 , 2009 ) . 
The authors show observational evidence of an emerging helical 
flux rope that apparently reconnects with the remnant magnetic 
field of a pre-existing filament. A Milne-Eddington (ME) inver- 
sion of the Stokes profiles observed with Hinode/SP over the re- 
gion where the transverse fields of the emerging flux rope were 
observed yielded a mean upflow of -0.3 ± 0.2 km s -1 . The ME 
inversion applied by these authors was able to separate the ve- 
locity of the magnetic and the non-magnetic components (see 
their Fig. 2). Interestingly, while the non-magnetic atmosphere 
showed a typical granulation pattern, upflows were systemati- 
cally inferred for the magnetic component. Thus, in that work, 
the picture of the filament formation is compatible with the exis- 
tence of transverse fie lds that are moving upward s into the chro- 
mosphere and corona. MacTaggart & Hoodl(l2010l) simulated the 
rise of a weakly twisted flux tube from the solar interior up into 



an overlying arcade in the corona. The axis of their flux rope 
became trapped slightly above the photosphere, while the field 
lines of the top of the flux rope were still able to reach coro- 
nal heights. In the same way, the 3D MH D simulation s of an 
emerging flux rope presented by lYelles Chaouche et al. 1 (120091) 
showed horizontal field lines that moved upwards with veloci- 
ties of ~ -0.67 km s~'. How ever, the autho r s poin ted out that 
the emergent flux studied by Oka moto et al.l (120091) was proba- 
bly less buoyant than the case simulated in their work, since the 
observed rising speed was smaller and the modification of the 
granular pattern le ss severe than in the simulations. Interestingly, 
Lites etal. (2010) also fo und a similar pattern to that presented 
bv lOkamoto et alJ y009) in their study of the evolution of a fil- 
ament channel inside an active region. The inferred magnetic 
structure in the filament channel corresponded to horizontal field 
lines elongated along the PIL (with equipartition strengths) that 
often (but not always) di splayed upflows (see the v mag panel in 
Fig. 5 of lLites et al.ll2010l) . A slow continuous upward motion in 
active region filaments had already been observed almost three 
decades ago ( Malherbe etal] 1 1 983t [Sciimieder et al.l 1 1 9851) . al- 
beit in chromospheric layers. It is thus clear that studying the 
behavior of transverse magnetic fields at the location of filaments 
(i.e., at AR PILs), is cr ucial to understanding th eir evolution (see 
also the discussion in IWelsch & Fisherl 12012b . In particular, it 
will help to shed some light on how the magnetic structures that 
sustain the filaments get filled with plasma and how they eventu- 
ally get rid of it, unloading most of their mass. To this end, it is 
crucial to obtain more data of the flows associated with filament 
formation and evolution. 

In Paper I, we presented the magnetic structure of an AR fil- 
ament, inferred from simultaneous spectropolarimetric ob serva- 
tions in the photosphere and chromosphere. The filament seemed 
to be divided into two parts: (1) one section that showed a flux 
rope configuration which had its main axis lying in the chro- 
mosphere; and (2) another section, which was likely to have a 
similar magnetic structure, but whose main axis was detected 
in the photosphere instead. This extremely low-lying part of the 
filament left an imprint on the photosphere in the shape of an or- 
phan penumbral system in the proximity of small pore associa- 
tions. Paper I suggested that these features were the photospheric 
counterpart of the filament. However, in order to understand the 
formation mechanism of this AR filament, a detailed study of 
the Doppler velocities, in a global scenario, needs to be carried 
out. The present paper addresses this issue by focusing in par- 
ticular on the photospheric velocities inferred from the Doppler 
shifts of the Si i line in order to search for systematic upflows and 
downflows. We also analyze the chromospheric Doppler shifts 
provided by the He i triplet. 



2. Observations 

The active region under study, NOAA AR 10781, was in its slow 
decay phase during our observing run in July 2005. The active 
region had a round leader sunspot of positive polarity and was 
followed by an extensive plage region that stretched towards 
higher latitudes. On top of the polarity inversion line (PIL), a 
filament with an inverse S-shape was unambiguously identified. 
Two spectropolarimetric maps were taken on July 3 (at latitude 
and longitude N16-E8, fi ~ 0.95), and seven on July 5 (at N16- 
W18, n ~ 0-91), usin g the Tenerife Infrared Polarimeter (TIP-II; 
Collad os et al . 2007) at the German Vacuum Tower Telescope 
(VTT) on Tenerife. In addition, a time series of ~ 19 minutes, 
with the slit fixed at the PIL, was acquired on July 5. 
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Fig. 1. Stokes / profile, normalized to the continuum intensity, 
observed by TIP-II near the polarity inversion line. The dashed 
vertical lines indicate the positions of the rest wavelengths for 
the spectral lines that are used in this work. 

The slit ((X'5 wide and 35" long, with a pixel size of (X'17 
along the slit) was aligned with the filament and a series of spec- 
tropolarimetric maps, centered on the filament, were acquired 
using a scanning step size of 074 or 073, depending on the see- 
ing conditions. Figure Q] shows the spectral range of the spec- 
tropolarimeter, which comprises the photospheric Si i 10827 A 
line, the chromospheric Hei 10830 A triplet and two telluric 
H2O lines. The helium triplet comprises a "blue" component and 
two blended "red" components. Table[T]provides the wavelength 
and the logarithm of the oscillator strength times the multiplic- 
ity of the lower level, log gf, for each line. This unique spec- 
tral window allows us to study the photosphere and the chromo- 
sphere simultaneously. The seeing conditions during the obser- 
vations were highl y variable and the adaptive optics system of 
the VTT (KAOS: fvon der Liihe etal]|2003h helped to substan- 
tially improve the image quality. The data sets were corrected 
for flat field, dark current an d polarimetric c alibration using the 
standard procedures for TIP (IColladosll 19991 120031) . 

Figure [2] shows an example, for both days, of the scanned 
area of the filament. From top to bottom, the slit-reconstructed 
images presented correspond to different wavelength bands: con- 
tinuum intensity, Sii line core intensity, Hei line core inten- 
sity centered at the red component and Si 1 line-of-sight mag- 
netogram. On July 3, the filament had a thin elongated shape 
in the chromosphere (see He 1 panel of Fig. [2), which we inter- 
pret to be the main axis or spine of the filament. Underneath, 
in the photosphere, the continuum image shows no outstanding 
features. The FOV of July 5 was slightly different and included 
newly appearing pores and orphan penumbral structures in the 
photosphere, mainly located along the PIL and below the fil- 
ament (see continuum image in the top righthand panel). The 
white dashed rectangles in Fig. [2] show the approximate portion 
of the filament common to both days of observation. The He 1 
intensity image of July 5 also shows the spine of the filament in 
the lower (common) part of the FOV. This feature is traced by 
an elongated shape in the Si 1 line core image, indicating that the 
structure lies at a lower mean height than on July 3. In the upper 
part of the He 1 image, on the other hand, the filament is more 
diffuse and less compact, and shows a series of dark threads that 
are highlighted by the white arrows (see also Fig. 5 of Paper I 
for more TIP-II images from July 5). Since the structure of the 
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Fig. 2. TIP-II slit-reconstructed images for both days of obser- 
vation. From top to bottom: continuum, Si i line core intensity, 
He i line core intensity (red component) and Si i Stokes V. The 
FOV was different on July 3 and July 5. We provide the dashed 
rectangles to show the approximate common area of the filament 
on both days. Characters a-d indicate the position of the Stokes 
profiles presented in Fig. [7] The white arrows point to high con- 
centrations of helium with elongated shapes (thread-like struc- 
tures). 



filament is identifiable in the silicon images, Paper I concluded 
that the filament was extremely low-lying and was directly re- 
sponsible for the orphan penumbral aggregations and pores. 
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Table 1. Atomic data used in this work. 



Line 


Wavelength (A) 


logg/ 


Si i 


10827.089" 


0.363* 


He i blue 


10829.091" 


-0.745+ 


He i red 


10830.250" 


-0.268 f 




10830.340" 


-0.047+ 


H 2 (telluric) 


10832.108* 





Notes. Wavelengths '"' from the NIST database {b] inferred in this work; 
see Appendix iBl 

References. f*1 lBorrero et alJ d2003l) : (f) VALD: lKupka et alj fl999!) 



Table 2. Examples of the line-of-sight to local solar reference 
frame transformation of the magnetic field inclination and az- 
imuth. The examples were averaged over 1 1 (PIL) and 8 (spine) 
pixels, in addition to the original binning of the data (see Sect. 

ED. 









LOS frame 




Solar frame 


Day 


Location 


Ion 


r(°) 


<M°) 




To O <f>e O 


July 3rd 


Spine 


He i 
Si i 


101.5 
84.2 


93.0 
63.6 


— » 
— » 


116.1 85.3 
98.1 54.6 


July 5th 


PIL 


He i 
Si i 


91.6 
88.9 


140.4 
114.6 


— > 
— » 


74.0 125.7 
81.2 100.0 



The main purpose of this paper is to identify what the plasma 
motions associated with the filament are, both in the chromo- 
sphere and the photosphere. In particular, we concentrate on 
whether the field lines are moving upwards, which would in- 
dicate an emerging process, or downwards, which would reveal 
a submergence phenomenon. The position on the solar disk of 
this active region is quite close to disk center on both days. This 
presents an advantage for the study of line-of-sight (LOS) ve- 
locities, which should not differ too much from the vertical ve- 
locities (measured in a local solar reference frame). For a better 
understanding of the inferred velocities we will distinguish be- 
tween two regions in the filament: the spine area and the diffused 
filament area (above the orphan penumbrae), the latter being ob- 
served only on July 5. As stated in Paper I, these two areas cor- 
respond to the same filament, but observed at different heights. 
It is therefore helpful to study them separately. 

A detailed analysis of the evolution of this AR was presented 
in Paper I. 

3. Data analysis 

To infer the physical parameters from the available sets of spec- 
tropolarimetric data we used two different inversion codes. It is 
important to recall that an inversion code provides full vector 
magnetic fields, but only the line-of-sight component of the ve- 
locity, in the observer's reference frame. Since one infers the 
three components of the vector magnetic field, the resulting in- 
clinations and azimuths can be projected onto the local solar ref- 
erence frame (i.e., with respect to the solar vertical, z, and along 
the solar latitude and longitude planes). However, difficulties 
arise when this transformation is carried out. The well-known 
180° ambiguity (due to the angular dependence of Stokes Q and 
U with the azimuth) provides two solutions for the magnetic 
field, with different azimuths and inclinations when projected 
onto the local solar reference frame. To find the correct solution, 
we used the AZAM code for t he disambiguation of the magnetic 
field vector dLites et al.l lY995) (see also Sect. 3.2 of Paper I for 
a more detailed explanation). Doppler shifts provide only one 
component of the velocity vector, the one projected along the 
LOS, rendering the projection onto the local solar frame impos- 
sible. Since the aim of this work is to compare the LOS velocities 
and their associated vector magnetic fields, for the remainder of 
this paper we will leave the magnetic field inclinations and az- 
imuths in the observer's frame. As our observations were done 
relatively close to disk center, the distinction between the two 
frames is more quantitative than qualitative. In order to show the 
impact of using the vector magnetic field in the observer's frame 
instead of in the local solar frame, Table [2] provides two exam- 
ples of the inclinations (y — > y ) and the azimuths (<p — > <p Q ), 
averaged over several pixels, in both reference frames. An incli- 



nation of y — 90° in the LOS frame means that the field lines 
are perpendicular to the LOS, while y Q = 90° means that they 
are parallel to the solar surface. As can be deduced from Table 
[2] while the exact values differ between the two frames, the hor- 
izontal or vertical nature of the field vectors is mostly preserved. 

3.1. Helium 10830 A and silicon 10827 A inversions 

Before inverting the helium spectra we carried out a 3 x6 binning 
of the data along the scanning and slit directions, respectively, 
and a binning of 3 pixels in the spectral domain to increase the 
signal-to-noise ratio (S/N), rendering a value of ~ 2000 for all 
maps. The resulting spectral sampling is ~ 33.1 mApx -1 and 
the pixel size 1.2 x 1.0 arcsec 2 and lxl arcsec 2 , for July 3 
and 5, respectively. However, since the Si i line is much stronger 
than the He i triplet, the inversions of the former were carried out 
without any binning, the pixel size being 0.40 x 0.17 arcsec 2 and 
0.30 x 0.17 arcsec 2 , on July 3 and 5 respectively. This preserved 
the original spatial resolution at an expense of a lower S/N of 
~ 500. The binned silicon data and inversions were used solely 
for the purpose of Table [2] 

The chromospheric Hei 10830 A triplet was inverted using 
a Mil ne-Eddington inversion code (MELANIE; Socas-Navarrol 
HEED). This code assumes that the source function varies linearly 
with optical depth while the remaining semi-empirical parame- 
ters stay constant with height. MELANIE computes the Zeeman- 
induced Stokes pr ofiles in the incomplete Paschen-Back (IPB) 
effect regime (see Socas-Na varro et al .1 120041 ISasso et al.l f2006. 
for a better understanding of the IPB effect on the Stokes pro- 
files in the Hei 10830A region) and has a set of nine free pa- 
rameters, which are iteratively modified by the code in order to 
obtain the best match between the synthetic and the observed 
Stokes data. The macroturbulence and the filling factor, or frac- 
tion of the magnetic component occupied inside each pixel, were 
fixed at 1.2kms _1 and / = 1 in the inversion, so no stray-light 
was used. Atomic-level polarization is not taken into account in 
MELANI E. However, this is not necessary for this AR filament, 
since [Kuckein et ail d2009l) proved that these He i Stokes profiles 
are dominated by the Zeeman effect. 

The photospheric Sii 10827 A line was inverted with the 
SIR code (Stokes Inversions base d on Response functions; 
iRuiz Cobo & del Toro Iniestalll992l) which solves the radiative 
transfer equation assuming local thermodynamic equilibrium 
(LTE). SIR gives a depth-dependent stratification of the inferred 
physical parameter as a function of the logarithm of the LOS 
continuum optical depth at 5000 A. 

We carried out two sets of inversions for the Si i line. The first 
set, referred to hereafter as "standard" inversions, used a uniform 
weighting scheme for the four Stokes parameters: wi^uy = 1- 
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An average intensity profile of the non-magnetic areas of each 
map (regions where Stokes Q, U, and V had negligible values), 
was used as representative of the stray-light. The stray-light frac- 
tion, a, used here as a proxy for the filling factor, was left as a 
free parameter in the inversion. Note that in the case of helium, 
the line is weak in quiet regions, so the use of a magnetic fill- 
ing factor does not provide any benefit. We studied the response 
function (RF) to velocity perturbations in several model atmo- 
spheres obtained from the inversion of the filament data, finding 
a consistent value of log t = —2 for the height at which the sen- 
sitivity was the highest. The height stratification derived by the 
inversion code was not used since the velocities at other optical 
depths had larger error bars. 

The Doppler shift of the stray-light profile cannot be set 
as a free parameter in the SIR inv ersion code. Thi s extra free- 
dom was used in the studies by lOkamoto et al.l d2009l) and 
iLites et all d2010l) to separate the magnetic and non-magnetic 
Doppler shifts. While the latter were dominated by granulation 
patterns, the former provided clear indications of upflows - a 
crucial result that we want to investigate further in this paper. To 
this end, we produced a second set of SIR inversions that were 
able to retrieve the Doppler shifts of the magnetic component 
without interference from a possible shift of the stray-light pro- 
file. These inversions, referred to heretofore as the "magnetic" 
inversions, forced a zero-weighting (wj = 0) for Stokes /, ba- 
sically eliminating any sensitivity to the intensity profile, and 
hence to the stray-light, in the inversion procedure. Furthermore, 
since we were interested in the Doppler shifts of the transverse 
fields near the PIL, we weighted Stokes Q and U equally with 
wq^u = 1 and set the weight of Stokes V to wy = 0.1. The filling 
factor was fixed to / — 1 . The remaining parameters were ini- 
tialized in the same way as in the "standard" silicon inversions. 

iBard & Carlssonl (120081) reported in their work that the Si i 
10827 A line is formed in non local thermodynamic equilibrium 
(NLTE) conditions. In Paper I we tested the NLTE effects on 
the inversions of the Si i line in the filament, using the departure 
coefficients, /3 (defined as the ratio between th e population den- 
sities c alculated in NLTE and LTE), reported in lBard & Carlssonl 
(2008) . We found almost negligible changes in the retrieved LOS 
velocities, temperature being the only parameter that was sub- 
stantially affected in the upper layers. Therefore, we decided 
not to consider the departure coefficients for the purpose of the 
present work. 

The reader is referred to Paper I for a more thorough descrip- 
tion of the inversions and NLTE effects. 

3.1.1. Two-component inversions 

The presence of several atypical Stokes V profiles (i.e., with sev- 
eral highly redshifted lobes) in the Si 1 10827 A line led us to use 
a two-magnetic component SIR inversion for these cases (see 
Sect. 14.21 ). Therefore, we had two initial guess atmospheres. For 
the first one, we used the same atmospheric model as that from 
the previous single-component inversion, but for the second, we 
increased the magnetic field strength from 500 to 1500 Gauss 
and the velocity from 0.5 to lOkms -1 . These initial conditions 
resulted in a good performance of the SIR code. 

3.2. Wavelength calibration 

The main aim of this work is to provide reliable Doppler veloc- 
ities measured at the filament at both heights. In order to cal- 
ibrate these velocities on an absolute scale, we had to correct 



them for several effects such as Earth's orbital motions and the 
solar gravitational redshift (see Appendix [A] for a complete de- 
scription of the calibration). The wavelengths presented in Table 
Q]were all obtained from the National Institute of Standards and 
Technology (NIST) except for the one corresponding to the tel- 
luric line. The literature reports several values, with differences 
of the order of tens of mA, which are large enough to cause 
systematic errors in our velocity calibration. We propose a new 
value for the rest wavelength of this telluric line, which we in- 
ferred by calibrating flat field TIP-II data taken at disk center 
(see AppendixIBl. 

4. Results 

4.1. Doppler velocities 

The inferred line-of-sight velocities of four representative cases 
(one on July 3 and three on July 5) are shown as slit- 
reconstructed maps in Fig. [3] Blue means upflow (v < 0) while 
red indicates downflow (v > 0). The pixels enclosed by the con- 
tours have LOS magnetic field inclinations between 75° < y < 
105°, y being the inclination inferred from the corresponding 
silicon or helium inversions. This criterion was chosen to study 
the velocities associated to transverse fields which, as reported 
in Paper I, dominate along the filament. The contours denote 
the area they belong to. The green contours enclose pixels with 
transverse field lines which lay at the orphan penumbral region 
whereas the black contours correspond to the spine region. Bear 
in mind that the FOV on July 5 only overlaps approximately with 
the upper half of the July 3 map (the white dashed rectangles in 
Fig.|2]give an idea of the common FOV). Moreover, small shifts 
of the FOV (in the vertical direction) between different maps of 
July 5 are also present. 

The panels in the upper row of Fig. [3]correspond to the Si i 
10827 A velocities, scaled between +1 kms _I , inferred from the 
standard inversions, i.e., those with Wi,q,u,v = 1. We found al- 
ternating blue and red patches in all panels, which revealed a 
mixture of upflows and downflows along the PIL, i.e., inside the 
contours. However, it is clear that on July 3, the inside of the 
black contour is more dominated by upflows than on July 5. On 
this day, there are no significant differences between the velocity 
patterns of the orphan penumbrae and the spine. 

We now turn to the panels of the middle row of Fig. [3] The 
velocity maps correspond to the magnetic inversions as defined 
in Sect. l3.1l (wf n nv = 0, 1 , 1,0.1). We therefore assume that the 
velocities in these panels are more likely to represent the motions 
of the transverse magnetic field lines (as opposed to the maps 
from the first row, which are more influenced by unmagnetized 
regions and flows along vertical field lines). The inversion code 
had problems in converging properly when Stokes Q and U were 
small, i.e., where the magnetic fields were weak or almost ver- 
tical. Nevertheless, along the PIL, the inversion code performed 
well and the uncertainties obtained were of the same order of 
magnitude as those from the standard inversions. On July 5, a 
comparison of the middle with the upper panels reveals that the 
"magnetic" Si i velocities inside the contours are globally more 
dominated by blueshifts, rather than redshifts. We interpret this 
result as an indication of a predominant upward motion of the 
transverse field lines in the photosphere, in particular inside the 
orphan penumbrae (green contours). Similarly, on July 3, there is 
more of an upward trend in the motions of the transverse fields. 

The velocity maps inferred from the (binned) He i 10830 A 
inversions are presented in the lower row of Fig. [3] Note that the 
color is now scaled between +2kms~ 1 . The first striking result 
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Jul. 3 14:39-15:01 Jul. 5 8:08-8:28 Jul. 5 9:33-9:52 Jul. 5 13:39-13:58 




10 20 30 5 10 15 20 35 5 10 15 20 25 5 10 15 20 25 

arcsec arcsec arcsec arcsec 



Fig. 3. Line-of-sight velocities inferred from: the Si 1 10827 A full Stokes "standard" inversions (top row), the Si 1 10827 A "magnetic" 
inversions (with weights: w/ = 0, Wqjj — 1 .0 and wy — 0. 1) (middle row) and He 1 10830A inversions (bottom row) at different times. 
Note that the color scale is saturated at [-1,1] and [-2,2] kms~' for silicon and helium respectively. Pixels inside the contours have 
magnetic field inclinations, with respect to the LOS, between 75° < y < 105°. In the July 5 maps, the colored contours distinguish 
between the orphan penumbrae area (green) and the spine region (black). On July 3, only the spine of the filament was seen. Negative 
velocities (blue) indicate upflows while positive velocities (red) indicate downflows. 



are the ubiquitous red areas (indicative of downflows) that ex- 
tend almost everywhere, with the exception of the region near 
the PIL. On July 3, an almost perfect correlation between the 
transverse field lines inferred from the helium inversions and the 
upflow areas can be seen within the black contour in the lower 
lefthand panel. Bearing in mind that on July 3 we clearly saw 
the filament axis in the helium core absorption map (see Fig. 0, 
and in Paper I we found sheared field lines parallel to its axis, 
one can easily deduce that the filament axis is rising in the chro- 
mosphere. Typical rising speeds are in the range [0,-1. 5]km s _1 . 

Two days later, on July 5, the spine region (black contours) 
hardly shows any upflows. The large-scale redshifted pattern 
seen elsewhere is, however, weaker. This means that, even if the 
filament axis was not moving upwards on July 5, it was still able 
to interfere with the mechanism that generates the large scale 
downflows (see Sect. 14.3b . In the diffuse filament region (above 
the orphan penumbrae, marked with green contours), the veloc- 
ity pattern is fundamentally different. Once again, the transverse 



fields near the PIL harbor clear signs of upflows. It is thus evi- 
dent that, while the general trend observed in He i corresponds 
to a large scale downflow, the filament axis (near the PIL) is a 
place where the downflows are not as strong and sporadic up- 
flow patches are also present. We will discuss this issue later, in 
Sectionl4~3l 

To understand the general trends in the velocity maps of Fig. 
[3] a statistical study of the average velocities (< v >) inside the 
black and green contours was carried out and compiled in Table 
[3] The number of averaged points (#) and the standard deviation 
(cr) are also provided in this table. We have found the following: 

1. On July 3, the mean velocities inferred from the Sii in- 
versions at the spine (i.e., below the filament axis) show 
an obvious blueshift. The magnetic inversions, averaging 
< v^ 38 >~ -0.150kms _I , are more blueshifted than the 
standard inversions. Upflows are also found in the chromo- 
sphere (< v He > 0.240 km s -1 ). 
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Table 3. Mean Doppler velocities (< v >), of all available maps for both days, calculated inside the contours, where the inferred 
LOS inclinations are between 75° < y < 105°, i.e., along the PIL. We distinguish between two different areas: the filament axis or 
spine (black contours in Fig. [3} and the orphan penumbral area (green contours), cr stands for the standard deviation and # is the 
number of pixels in each area selected for the study. 



Spine Orphan penumbrae 



Day 


Time range 


Inversion 


< v > 


cr 


# 


< v > 


(X 


# 




(UT) 




(kms- 1 ) 


(kms" 1 ) 


points 


(kms" 1 ) 


(kms- 1 ) 


points 






Si 


-0.123 


0.725 


2097 








3 Jul. 


13:53 - 14:14 


Si ma 8(t) 


-0.134 


0.546 


1974 


No orph 


an penumb 


'ac 






He 


- 0.266 


0.574 


146 












Si 


-0.010 


0.542 


1114 








3 Jul. 


14:39-15:01 


Si raag 


-0.163 


0.406 


1088 


No orph 


an penumb 


'ac 






He 


-0.216 


0.521 


94 












Si 


- 0.011 


0.424 


1019 


+0.086 


0.383 


831 


5 Jul. 


7:36-8:05 


gjmag 


- 0.190 


0.333 


1018 


-0.250 


0.329 


825 






He 


+ 0.808 


0.566 


42 


+0.547 


1.303 


54 






Si 


+ 0.108 


0.326 


870 


+0.026 


0.353 


1009 


5 Jul. 


8:08 - 8:28 


Si raag 


- 0.144 


0.263 


870 


-0.241 


0.368 


997 






He 


+ 0.884 


0.604 


42 


+0.346 


0.951 


59 






Si 


- 0.011 


0.186 


780 


+0.079 


0.317 


1127 


5 Jul. 


8:42-9:01 


Si mag 


- 0.237 


0.288 


780 


-0.243 


0.353 


1106 






He 


+ 0.954 


0.564 


36 


+0.427 


1.002 


63 






Si 


- 0.040 


0.267 


760 


+0.110 


0.322 


1023 


5 Jul. 


9:02-9:21 


Si mag 


- 0.291 


0.246 


760 


-0.229 


0.336 


1006 






He 


+ 0.649 


0.464 


37 


+0.183 


0.622 


62 






Si 


+ 0.045 


0.263 


569 


+0.021 


0.332 


1497 


5 Jul. 


9:33-9:52 


Si mag 


- 0.085 


0.248 


569 


-0.209 


0.400 


1457 






He 


+ 1.013 


0.645 


22 


+0.366 


0.826 


102 






Si 


-0.127 


0.282 


755 


+0.021 


0.314 


1292 


5 Jul. 


13:39-13:58 


Si mag 


- 0.063 


0.248 


755 


-0.283 


0.419 


1273 






He 


+ 0.779 


0.960 


41 


+0.278 


0.814 


60 






Si 


- 0.046 


0.375 


618 


+0.211 


0.301 


1257 


5 Jul. 


14:31 -14:51 


Si mag 


+ 0.014 


0.290 


618 


-0.099 


0.427 


1248 






He 


+ 1.083 


1.009 


30 


+0.271 


0.558 


75 



Notes. (l) Velocities above 2 kms 1 or below -2 kms 1 inferred from the "magnetic" Sii inversions (and mostly due to weak Q and U signals) 
were excluded from the statistics. However, only very few points were affected by this filter ( between 9 - 40; and 123 for the first map). 



2. For the first four maps of July 5, the spine presents an upward 
motion of the photospheric transverse fields (black contours) 
as inferred from the magnetic inversions, with typical veloc- 
ities of < v™ s >= -0.215kms _1 . In the last three maps 
the velocities drop to almost zero. The standard Si i inver- 
sions show, by contrast, velocities which oscillate around 
zero (< v si > 0.01210ns- 1 ). 

3. On July 5, the chromospheric velocities at the spine 
are dominated by downflows in the range of v'H e e 
[0.65, 1.08] kms -1 . We note, however, that these downflows 
are smaller than the typical values observed elsewhere in the 
surrounding facular region (see Sect. 14.3) . 

4. At the orphan penumbrae (green contours), there is a sys- 
tematic difference between the photospheric LOS veloci- 
ties inferred from the "standard" and those obtained from 
the "magnetic" inversions of the Sii line. While the lat- 
ter clearly show upflows in all maps, with values be- 
tween v™ dS e [-0.099, -0.283] km s _1 and an average of 

< v™ ; ag >= -0.222km s _1 , the standard inversions show 
slightly slower velocities, with downward motions in the 
range of Vsi e [0.021, 0.211] km s , and a mean value of 

< v si >= 0.079 km s -1 . 

5. Above the orphan penumbrae (inside the green contours) 
the average velocities inferred from the He i inversions show 
downflows in the range of vn e e [0.18, 0.55] kms -1 . These 



downflows are weaker than those observed in the spine (see 
point 3 above). 

4.2. Two-component Stokes V profiles at the PIL: supersonic 
downflows 

On July 5, the day when we first observed the orphan penumbral 
system, we detected the presence of atypical multilobed Stokes 
V profiles (e.g., Fig. |4|i in the Sii 10827 A line. These multi- 
lobed profiles are a clear indication of the presence of several 
magnetized components within the resolution element. Most of 
these multi-component profiles were found near the borders of 
the orphan penumbrae. To illustrate this, we selected one map 
from July 5 (taken between 13:39-13:58 UT) and marked the 
location of all two-component profiles with small cyan crosses 
(see Fig. [5]). The other maps from July 5 show a similar spatial 
distribution of the multi-component profiles. Only points with a 
Stokes V/I c signal of at least 0.01 were selected for this study. 
In this way, we made sure that the origin of these profiles was 
related to a multiple-velocity component scenario, and we were 
able to rule out other causes such as mixed polarities, which are 
also a very common occurrence in the proximity of PILs. Figure 
13 a) shows a continuum intensity image. The two-component 
profiles are located on top of pores and orphan penumbrae, be- 
tween x e [16", 22"]. The helium red component intensity map 
in Fig. |5jb) shows that the largest group of atypical Stokes V 
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Fig. 4. From left to right and top to bottom: Stokes /, Q, U and V profiles (dots) of the Si i line normalized to the continuum, and 
best fits (solid line) using a two-component SIR inversion. The profiles correspond to [x,y] = [20'.'1,19'.'4] (area 4 in Fig. [5}. The 
single-component fit (i.e., standard fit) is shown in the dashed lines for comparison and clearly does a worse job at fitting Stokes V. 
The physical parameters inferred from the two-component inversion are: Vi = 0.3 + 0.2 km s _1 , B\ = 1001 + 68 G, j\ = 67° + 3°, 
0i = 114° + 2°,/ ~ 0.94 for the standard component and v 2 = 11 ±2kms- 1 , B 2 = 1424 + 300G, y 2 = 23° +20°, 4> 2 = 81° +51°, 
/ ~ 0.06 for the magnetic one. The inferred stray-light is a ~ 5%. The uncertainties associated to the parameters of the second 
component were obtained with SIR and verified with the sensitivity analysis explained at the end of Section l4~2l 



profiles is co-spatial with the dark He i thread located around 
[x,y] = [20", 20"]. Only three asymmetric profiles were de- 
tected at the spine, around [x,y] = [18"5,3"]. The contours in 
Figs-Hb-c) correspond to the darkest areas, i.e., highest absorp- 
tion of the He i red component, of Fig. [2b). Bear in mind that the 
multi-component profiles were only found in the photosphere; 
there is no signature of these highly Doppler shifted components 
in the chromospheric He 1 10830 A profiles. The atypical Si i pro- 
files contain information about a strongly redshifted component 
that the standard SIR inversions cannot pick out. The conspic- 
uous signature present in Stokes V calls for a two -atmosphere 
inversion that can provide some insights into the nature of this 
shifted component. Such inversions were performed on those 
specific profiles where this signature was visually detected. The 
following properties were inferred: 

1. The atypical Stokes V profiles appear at, or very close to, 
the PIL, i.e., where the magnetic field lines are mainly hor- 
izontal (see inclinations inferred from the standard single- 
component silicon inversions, ysi in Fig. [3d)). 

2. From the two-component SIR inversions we inferred: (1) a 
dominant transverse component with a filling factor in the 
range of f\ e [0.87, 0.98] and inclinations with respect to 
the LOS n e [68°, 95°]. This is the component that is most 
similar to that retrieved by the standard inversions. (2) A sec- 
ond component with smaller filling factors, fa e [0.02, 0.13], 
oriented more parallel to the LOS, with inclinations y 2 e 
[21°, 58°]. 

3. The inferred Doppler maps show that the dominant trans- 
verse component has velocities in the range of v; e 
[-0.1,0.5] kms -1 while the longitudinal one, with smaller 



Table 4. Mean LOS inclinations (< y >) and velocities (< v >) 
inferred from the two-component inversions of the Si i line. The 
different areas are marked in Fig. HJc). cr stands for the standard 
deviation and # is the number of pixels averaged in each area. 



Area 


Comp. 


< y > 


<r y 


< v > 


<r r 


# 






o 


O 


(kms" 1 ) 


(kms" 1 ) 


(px) 


1 


1 


94.6 


1.4 


-0.126 


0.135 


25 


2 


57.7 


51.5 


11.890 


8.867 




1 


76.3 


2.3 


0.351 


0.034 




2 


2 


20.7 


20.2 


4.699 


2.439 


3 


3 


1 


68.5 


6.5 


0.496 


0.315 


32 




2 


43.7 


18.1 


7.595 


2.278 




1 


72.7 


11.8 


0.355 


0.256 


85 


4 


2 


28.4 


17.4 


9.632 


2.795 




1 


78.7 


3.5 


0.250 


0.113 


18 


5 


2 


49.4 


20.4 


6.190 


2.919 



filling factors, has downward velocities v 2 e [6.2,11.9] 
kms 1 . It is clearly a supersonic downflowing component. 
4. In Fig. [5j all the detected atypical Stokes V profiles had a 
second component with the same polarity than the first one, 
except for the group of points between y = [11", 13"] (area 
1 in Fig. [3c)) which had a second component with opposite 
polarity to that of the first one. Note, however, that the domi- 
nant polarity, being highly transverse, is very much prone to 
polarity changes due to projection effects. 

The multiple component profiles were grouped in five sepa- 
rate areas. Since the properties of the dominant and of the red- 
shifted components differ in each of these areas, a separate sta- 
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Fig. 5. (a) The continuum intensity image show pores and orphan penumbrae at the PIL. (b) The absorption in helium (red core) 
clearly shows the spine in the lower part of the image and a diffuse filament in the upper part, (c) LOS velocity map inferred 
from the standard single-component inversions of SIR. The marked areas (1-5) correspond to the two-component Stokes V profiles 
(the average inclinations and velocities for these areas are presented in Table 3J. (d) LOS inclinations from the single-component 
inversions. In all maps, the cyan colored crosses mark the positions where two-component profiles were detected. 




-1.0 -0.5 0.0 
AA [A] 



-1.0 -0.5 0.0 0.5 
AA [A] 




-1.0 -0.5 0.0 0.5 
AA [A] 



Fig. 6. Dots represent the observed Sii Stokes U, Q, and V profiles (same as Fig. |4}. The overplotted synthetic profiles were 
computed using different magnetic field inclination angles for the second component, that varied between y e [40°, 120°]. The 
test was performed by varying the magnetic field inclination of the second component and keeping the remaining atmospheric 
parameters (which came from the best fit of the SIR inversion) fixed. The dotted red line shows the best matching profile, whose 
corresponding atmospheric parameters are shown in the caption of Fig. [4] 



tistical study for them has been carried out. The five areas are 
marked in Fig. [3c ). Table [4] summarizes the mean LOS incli- 
nations (< y >) and velocities (< v >), for both components, 
in the five areas. The standard deviation (<x) of both quantities, 
as well as the number of pixels used for the statistics (#), are 
also listed. As stated before, the second component is completely 
dominated by magnetic field line inclinations which are more 
longitudinally orientated. These inclinations, however, vary sub- 
stantially within each area (note the large values of the standard 
deviation <x r for Comp. 2 in Table @}. In these atypical cases, 
the inversion code provides considerable uncertainties associ- 
ated with the atmospheric parameters, something that should be 
expected since this redshifted component has a small effect on 
the observed profiles, as reflected by its low filling factor (mostly 



< 0.1). However, a clear tendency for smaller inclinations, i.e., 
more vertically oriented field lines, is seen in all five areas. 

Given the small amplitude of the signatures that this com- 
ponent leaves on the emergent Stokes profiles, the question is 
raised of how reliable the inversions are. To address this ques- 
tion, we performed the following test: We selected several pix- 
els with multilobed Stokes V signals and used the atmospheres 
resulting from the two-component inversions to synthesize sets 
of Stokes profiles. For each model atmosphere, several different 
realizations of the synthesis, with different values of the inclina- 
tions and azimuths of the second component, were carried out 
(the first component was kept unchanged). Figure [6] shows one 
example of the synthetic Stokes Q, U, and V profiles (Stokes / 
showed no changes and was therefore not presented) with differ- 
ent inclinations, between 40° and 120°. Stokes Q and U present 
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Fig. 7. Left: averaged intensity profiles 7// c of four selected positions of the map acquired on July 5 between 13:39-13:58 UT. Each 
profile is marked with a character (a-d) which indicates its position in Fig. [2] Note the shifted red wing of the blended He i red 
component for the two facular profiles indicating the presence of more than one component contributing to the line formation. 
Right: the corresponding helium Stokes V/I c profiles to the intensity profiles shown on the left panel. The solid lines along the 
jc-axis indicate V/I c = for each profile. The dotted vertical lines, in both panels, mark the position of the line cores of the involved 
ions. 



minor changes in their shape that are just above the noise level 
and therefore not conclusive. However, the changes in inclina- 
tion substantially affected the resulting Stokes V profiles. It is 
apparent from Fig. [6] that the most likely inclination is far from 
being transverse. In fact, the best fit (dotted red line) yields an 
inclination of y — 22.9° + 20°. We also carried out the same 
test changing the azimuth. Stokes Q and U were only slightly 
affected by these variations, again within the noise level, while 
Stokes V was completely insensitive to the changes. This shows 
that while the inclination of the redshifted component is well 
established (and that it is more vertical than the dominant, trans- 
verse, component), the value of its azimuth is highly uncertain. 



4.3. Ubiquitous downflows around the PIL in the helium 
velocity maps 

As mentioned above, the velocities inferred from the He i inver- 
sions outside the PIL (in the faculae) are ubiquitously dominated 
by downflows on both days (see the lower panels of Fig. [3}. This 
was largely unexpected, so it is important to understand what 
signature in the spectral profiles made MELANIE infer these 
downflows. Figure [7] shows four representative Stokes I and V 
profiles from different areas of the map taken between 13:39- 
13:58 UT on July 5. They correspond to: (a) the spine, (b) the 
PIL, (c) the negative and (d) the positive polarities within the fac- 
ulae. The locations of pixels a-d are marked in Fig. [2] The dotted 
vertical lines, in Fig. [7] mark the centers of the line cores at rest. 
In the lefthand panel, the intensity spectra comprise the photo- 
spheric Si i line and the He i triplet. The cores of the four silicon 
intensity profiles are almost perfectly centered at their rest wave- 
lengths. However, while the Hei triplet profiles of the spine and 
the PIL are also approximately at rest, the two facular profiles are 
clearly redshifted. It is also of interest that the red wing of the 
He i red component of the facular profiles is stretched away from 
the line core. This, together with the fact that the profiles are not 
symmetric, indicates that more than one atmospheric component 
was involved in the line formation process, and that one of these 



components harbors distinct downflows. The redshift in the he- 
lium red component of the facular regions is also seen in their 
Stokes V profiles (c and d) in the righthand panel of Fig. |7j As 
expected, profile (b) has almost no Stokes V signal since it is 
located at the PIL, where the transverse fields dominate. 

It is, thus, clear that the facular profiles seen in the chromo- 
sphere are substantially redshifted. A possible explanation for 
this is given in Sect. [5] To quantify the general trend of the facu- 
lar velocities found in our data sets we did averages only taking 
into account those pixels which had photospheric inclinations 
close to longitudinal (y Si < 25° and y S i > 155°). Note that the 
inclinations we used are the ones inferred from the Si i 10827 A 
inversions, whereas the velocities correspond to the ones ob- 
tained from the Hei 10830 A inversions. This criterion ideally 
distinguishes the points belonging to the fac ulae (having mainly 
longi tudinal fields at photospheric heights, Mart inez Pillet et aTl 
[199% from those of other areas, such as the PIL or outside the 
faculae. The mean velocities found in the faculae, for all maps 
(July 3 and 5), were in the range of vgf ~ 1.4 - 1.8 kms -1 . The 
standard deviations associated to these velocities were typically 
around cr - 1.0 - 1.4kms _1 , and the number of points used for 
the statistics of each map was between 300 - 400. As mentioned 
above, since these He i intensity profiles are likely to comprise 
more than one comp onent, multi-component inversions in this 
spectral region (e.g., lLagg et aDl2007t ISasso et al.ll201 ll) might 
reveal slightly stronger downflows in the faculae. 

4.4. Time series 

A time series with the slit fixed over the PIL was acquired on 
July 5, approximately between x ~ 18" - 20" on the righthand 
panel of Fig. [2] The velocities shown in Fig. [8] were retrieved 
from the standard Si i inversions. Time is represented along the 
x-axis while the y-axis shows the spatial direction along the slit. 
In the lower half of the figure, the photospheric 5 minute oscilla- 
tion is unequivocal. This corresponds to the area below the spine. 
There is an apparent phase shift of the oscillation pattern, start- 
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ing at y ~ 15". This shift coincides with the beginning of the 
orphan penumbrae. The single most striking finding from this 
figure was the continuous downfiow seen between y ~ 23" -25" 
(dark red pattern) throughout the whole time range, with which 
not even the photospheric oscillation pattern interfered. 

The Stokes V profiles were carefully inspected one by one 
along the slit every minute. This was done in order to find the 
atypical multilobed Stokes V profiles defined in Sect. 14.21 As in 
Fig-El the cyan crosses in Fig.|8]represent the location of the de- 
tected atypical Stokes V profiles. From this figure it is obvious 
that the lifetime of the supersonic downfiow was at least as long 
as our time series (19 minutes). We take this time to be a lower 
limit but it might, in fact, be much longer, since this supersonic 
downfiow was not an isolated event. Other strong downflows 
have been detected on all maps from July 5 (albeit not contin- 
uously), starting from 7:36 until 14:51 UT. 



o 20 



td 15 




10 

Time [min] 

Fig. 8. Photospheric velocity time evolution. The time series was 
taken with the slit fixed on top of the PIL (between x ~ 18" -20" 
in the second column of Fig. [2] ). The velocities were inferred 
from the "standard" Sii 10827 A inversions. The cyan-colored 
crosses mark the positions of the two-component profiles. 



5. Discussion 

In this paper, we have investigated the line-of-sight velocities in- 
ferred from chromospheric and photospheric spectropolari met- 
ric inversions of the Stokes profiles measured in and below an 
active region filament. The aim of this study, together with the 
results of Paper I, is to clarify the formation and evolution of 
this AR filament. In Paper I we reported that the magnetic field 
configuration of the spine of the filament is compatible with a 
flux rope topology (on July 3, and also in the spine portion of 
July 5). Furthermore, the orphan penumbral region of July 5, 
was also likely to be the imprint of a flux rope, but with its 
axis sitting low in the photosphere. These results are supported 
by recent magnetic fie l d extr apolations of the same filament by 
Yelles Chaouche et al. (2012) who, for July 5, also deduced the 
presence of a flux rope structure whose axis lay below the for- 
mation height of the Hei 10830 A lines. Our AR filament is, 
therefore, extremely low-lying. 

The next question that needs to be addressed is the mecha- 
nism by which flux ropes in AR filaments are formed. First of all, 
at the time of the observing run, the filament had already reached 



the chromosphere. However, on July 5, the low-lying magnetic 
structure and the newly appeared pores and orphan penumbrae 
indicated that new flux had emerged between the two days of 
observation. This flux was mainly horizontal due to its proxim- 
ity to the polarity inversion line, as reported in Paper I. Thus, 
we placed special emphasis on the study of the photospheric 
LOS velocities associated with these transverse fields, i.e., in- 
ferred from what we called the silicon "magnetic" inversions, 
which provided information similar to th at of the Dopplergrams 
of the ma gnetic comp onent presented by Okam oto et al.1 (2009) 
and iLites et al.l d2010h . There is a very good agreement among 
the velocities from all maps of July 5 (see Table [3). All the 
maps but one pointed to an upward movement of the pho- 
tospheric transverse fields, with average velocities between < 
v™ ag >= -0.21 and -0.28 km s~', the last map being the excep- 
tion, where slower upward motions (~ -O.lOkms -1 ) are shown. 
Th ese velocities agree su rprisingly well with the ones presented 
bv lOkamoto et al.l l2009) (-0.30 + 0.20km s _1 ). However, when 
carrying out "standard" inversions, (weighting all four Stokes 
parameters in the inversion equally) the average inferred veloci- 
ties were close to zero or positive (downflows), with values be- 
tween < vsi >= 0.02 - O.llkms -1 and 0.21kms _1 for the last 
map. Note that the differences between both inversions is, on av- 
erage, of around 0.3kms~'. This is indeed a small number. It 
corresponds to 11 mA or, equivalently, a sixth of the spectral 
resolution. The fact that we obtain velocities that are more red- 
shifted when we weight Stokes / and V equally in the inversion 
indicates that a longitudinal component harboring downflows is 
present, to some degree, in all pixels. While no evidence for this 
is provided here, we hypothesize that the same chromospheric 
downfiow that is observed almost everywhere in He i has an im- 
pact on the Si i data, producing slightly more redshifted values 
when all Stokes parameters are equally weighted in the inversion 
code. 

In the chromosphere above the orphan penumbrae, the re- 
trieved mean LOS velocities are dominated by downflows in the 
range of < ve e >= 0.18 - 0.55kms~'. Nevertheless, these av- 
eraged velocities need to be interpreted with caution since the 
helium velocity maps of Fig. [3] clearly show localized upflows 
and slower downward motions at the PIL (inside the green con- 
tours). These upward motions at the PIL reach velocities of about 
-1 km s in all the maps of July 5. A continuous upflow of mass, 
lasting at least ~ 7 hours, is detected at some locations of the 
PIL. Since helium traces the upper part of the flux rope (whose 
axis is lying below, in the photosphere), these localized upflows 
strongly suggest that the material is being pushed upwards by 
the magnetic field lines of the to p of the flux rope (similar t o 
what is seen in the simulations by Mart fnez-Svkora et al.l l2008). 
This is in ag reement wi t h previous simulations of flux rope emer- 
gence (e.g.. lFan|l200lt Manchester et al.ll2004t lArchontis et al.1 
2004; Fan 2009), where only the upper field lines are able to 
rise and expand in the corona, leaving the main body of the flux 
tube behind. This process is driven by the so-called Parker insta- 
bility. This interpretation would lead to plasma drainage wher- 
ever the field lines are not horizontal or dipped, and, hence, 
would explain the photospheric downflows close to the PIL. 
Sometimes, the downflows at the orphan penumbrae produce 
atypical multilobed Stokes V profiles. Examples of these were 
found in some areas near the PIL (see Fig. |5j- Two-component 
inversions of these Stokes V profiles yielded mainly supersonic 
speeds, in the range of 6 - 12kms _1 , for the redshifted compo- 
nent. Similar strong downflows have been detected in the past 
re.g JMartinez Pillet et al.l 1 994t Otes et al.l2002tlShimizu et all 
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2008) and have al so been found in nu merical models of emerg- 
ing flux regions (ICheung et al.l 1200 81. The case described by 
Marti nez Pillet et al.l d 19941) is particularly relevant to this paper 
because it describes a downflow below an AR filament, similar 
to the ones studied here. Magnetic reconnection is one of the 
proposed mechanisms to cause th ese supersonic down flows as 
suggested from the simulations of Cheung et al. (2008). In our 
data sets, this could only have taken place in the photosphere 
which is where we see them clearly. The long periods of time 
(tens of minutes) during which the redshifted flows persist ba- 
sically unaltered, almost rule out any link with a fast process 
such as reconnection. The downflows might also be interpreted 
as a combination of plasma draining from the flux rope's body 
during its emergence process at the photosphere and/or the con- 
tinuation of the ubiquitous chromospheric downflows. 

The filament forma tion model of 

van Ballegooii en & Martens! dl989l) would expect hori- 
zontal field lines s ubmerging at the PIL after a magnetic 
reconnection event (van Ballegooiien & Mackav 2007). While 
strong downflows in our observations are clearly detected at 
several positions along the PIL, they are never co-spatial with 
transverse magnetic fields. The inclinations retrieved from the 
inversions of the atypical Stokes V profiles, showed that the 
second component (the strongly redshifted one) was much 
more longitudinally orientated (72 e [21°, 58°]) than the first 
one (71 6 [68°, 95°]) (see more details in Table g). Our test 
inversions support the reliability of these results (see Fig. [6]). 
Therefore, we can safely say that our AR filament presents no 
evidence of submerging horizontal field lines, thus, challenging 
some features of the aforementioned models. 

Let us consider now the spine of the filament, which was 
seen in the middle of the FOV on July 3 and only in its lower 
part on July 5. This portion of the filament is what we would call 
a classical case of AR filament, "classical" meaning that it has a 
thin elongated shape presenting a flux rope topology and is lo- 
cated in the chromosphere (it is perfectly discernible in the He 1 
absorption image of Fig. There is a distinct difference be- 
tween the chromospheric LOS velocities of both days. On July 
3, the spine is completely dominated by upflows which perfectly 
trace the transverse fields (see black contour in the lower left 
panel of Fig. |3J. The average upward velocity inside this black 
contour was < v'H e >~ -0.24 km s -1 . On July 5, however, the 
spine is not as strongly blueshifted, characteristic that contin- 
ues throughout the day (see Fig. [3}. Downflows in the range of 
< >~ 0.81 and 1.09kms _I , are detected in the ~ 7 hours 
between the first and last maps. We interpret this result as a halt 
in the rise of the filament's axis. This event makes the spine area 
be dominated by the same ubiquitous chromospheric downflows 
that are seen in the rest of the FOV. 

The photospheric motions of the spine on July 3 show the 
same trend as the chromospheric ones, with average upflows of 
the transverse fields of < v'™ g >~ -0.15kms _1 . We interpret 
this as a coordinated upflow of the entire flux rope. This is com- 
patible with the flux rope emergence simulations, in which the 
whole structure is very low-lying. On July 5, the transverse fields 
begin the day moving upwards < v™ 8 >= -0.19kms _1 , but af- 
ter ~ 2 hours they seem to stop, still remaining at rest by the end 
of the observing run, four hours later. 

We found ubiquitous downflows in the chromosphere, on 
both sides of the polarity inversion line. This was consis- 
tent throughout all of our observations. The average line-of- 
sight velocity in the faculae, for all maps, was 1.6kms 
with a dispersion of ~ 1.2kms _1 . We hypothesize that 
these downflows are a manifestation of the so-called coro- 



nal rain. This phenomenon occurs when plasma condenses 
in the corona and then flows along coro nal loops, down 
into active regions ( Tandber g-Hanssenll 1 9951) . So far, however, 
our knowledge of this process is sti ll rather poor (but see 
lAntolin & Rouppe van der Voortl I20I2I Recent studies reveal 
that the velociti es found in coronal rain are in a range between 
20- 120 km s" 1 (lAntolin et alJ2010t lAntolin & Verwichtd201 it 
lAntolin & Rouppe van der Voortl2 012). This is much faster than 
the velocities obtained from our data (Fig. [3}, although our mea- 
surements correspond to lower heights. There is a strong pos- 
sibility that more than one atmospheric component, along the 
line-of-sight, contributes to the Hei 10830 A triplet formation. 
This can be deduced from the shape of the red wings of the in- 
tensity profiles in the faculae (lefthand panel; Fig . [7j) . However, 
the signature is weak and, hence, it is unclear whether strong 
downflows, such as those re ported by previous au t hors in the He 1 
triplet (~ 42kms~ 1 , e .g., Muglach etail 119971: ISchmidt etakl 
200d: lLagg et~alll2007l) . would be inferred in a multi-component 
analysis. Another explanation could be that material is falling 
from the slowly rising flux rope structure on July 3, or from 
the expanding field lines above the orphan penumbrae on July 

5, following vertical magnetic field lines (in a similar way as the 
draining of rising loops proposed by Lagg et al. 2007). However, 
these processes can also be considered to be different manifesta- 
tions of the coronal-rain phenomenon. 

6. Conclusions 

For the sake of clarity we divided the study of the AR fila- 
ment into two parts, the main distinction between them being 
the height at which the filament axis lies. The first part corre- 
sponds to the spine (seen mainly on July 3, and also in the lower 
half of the FOV on July 5), where the filament axis lies in the 
chromosphere. The second corresponds to the diffuse filament 
(seen only in the upper half of the FOV on July 5), which sits 
above the orphan penumbrae and pore regions, and has a much 
lower-lying filament axis. In this latter case, the helium 10830 A 
only traces the upper part of the flux rope, as explained in Paper 
I. The main conclusions of this study are: 

1. On July 3, the LOS velocities inferred from the helium and 
"magnetic" silicon inversions (that trace the behavior of the 
horizontal magnetic fields) in the spine region, show general- 
ized upflows, which we interpret to represent the emergence 
of the flux rope structure as a whole. Two days later, on July 
5, the spine shows downflows in the chromosphere, similar 
to those seen elsewhere in the facular region. Yet the photo- 
spheric velocities in this region present upflows for the first 
few hours of the day that drop to velocities close to zero to- 
wards the end of the observing run. 

2. On July 5, in the chromosphere, above the orphan penum- 
brae, the LOS velocities are on average redshifted. This red- 
shift is, however, smaller than the dominant redshift seen in 
the facular region. Indeed, blueshifted patches are present 
along the PIL in all data sets. We propose that the blueshifted 
patches seen in the chromosphere are due to field lines which 
expand from the lower-lying flux rope into the chromo- 
sphere, similar to what is found in the flux rop e emergence 
simulations (e.g.. iMartfnez-Svkora et al1l2008l) . The photo- 
spheric transverse field lines along the orphan penumbrae are 
clearly moving upwards. This is consistent in all seven data 
sets. Therefore the sheared field lines (the flux rope axis) are 
rising. 
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3. Atypical multilobed Stokes V profiles were found in the 
photosphere near the PIL. Two-component inversions of 
these profiles revealed localized supersonic downflows in 
the strongly redshifted component. These downflows last 
for at least 19 minutes, ruling out any episodic origin. 
Furthermore, the retrieved magnetic fields associated with 
this component are oriented along (or very close to) the 
line of sight. Therefore, we cannot identify this compo- 
nent with submergi ng loops that harbor horizon t al fie lds as 
in the models fromlvan B allegooii en & Martensl (1 19891) and 
Ivan B allegooiien & Mackav (2007]). It is also important to 
point out that the number of detected two-component profiles 
was only a small percentage of the total number of points in- 
side the PIL (12 %), which, on average, showed an upward 
motion of the transverse fields. 

4. Almost ubiquitous redshifts of the chromospheric Hei 
10830 A triplet, with average downflows of ~ 1 .6 km s^ 1 and 
a dispersion of ~ 1.2kms _1 , were found in the faculae in all 
data sets. 

The global picture resulting from the Doppler shifts stud- 
ied in this paper can be summarized as follows. AR faculae are 
immersed in a global rainfall of mass from the upper layers, as 
indicated by the redshifted He i line. We have tentatively associ- 
ated this redshift with coronal rain. This process has no evident 
counterpart in the photosphere, being some localized supersonic 
downflows the only possible candidate. In this environment, the 
transverse field lines (including the filament axis observed ei- 
ther in the chromosphere or in the photosphere) display upflows. 
The only case where this does not hold true is that of the chro- 
mospheric spine region on July 5. In this area, only the global 
downflow is observed. 

As we have learned from the simulations, the process of flux 
rope emergence through the photosphere is not an easy one and 
will often be aborted. We propose that this is what might have 
happened on July 5 at the spine region. Note, however, that the 
rest of the filament for this day displayed clear upflows at both 
heights. All in all, the present study supports the scenario of an 
emerging flux rope from below the photosphere. Vector magne- 
tograms, as well as LOS velocities in the photosphere and chro- 
mosphere, agree with the proposed scenario. Whether or not this 
emergence process is common to all active region filaments, or 
at least to those that present a photospheric manifestation (or- 
phan penumbrae), needs to be proven by using further multi- 
wavelength and multiheight observations. 

There are some limitations in this work that need to be con- 
sidered. First, an important observational gap on July 4 exists. 
Data from that day would have contributed to a better under- 
standing of the inferred magnetic structure observed on July 5. 
Another limitation to this work is the small field of view of the 
polarimeter that we had at the time of the observations (nowa- 
days the slit is twice as large). This work has proved that there is 
an imperative need for multi-wavelength instruments and larger 
FOVs in order to understand the formation process and evolution 
of AR filaments. It is crucial to have magnetic field information 
of at least two heights, e.g., in the photosphere and in the chro- 
mosphere, to carry out a simultaneous and co-spatial analysis of 
the evolution of the filament. 

The flux rope that constitutes the studied AR filament is ex- 
tremely low-lying, especially the portion observed on the sec- 
ond day. This result makes it hard for the filament formation 
models tha t build fl ux ropes in the corona by reconnection (e.g. , 
Ivan Ballegooiien & Martensl 1 1989t iDeVore & Antiochosl 12000) 
to reproduce the scenario that we observe. In particular, the large 



scale submergence of photospheric transverse field lines is not 
observed in our data. While a rather complete picture of the evo- 
lution of this AR filament (favoring a flux rope emergence from 
below the photosphere) has come out of this series of papers, we 
would like to stress that the formation process of AR filaments 
in general could differ, maybe substantially, from this case. 

Appendix A: Velocity calibration 

In order to obtain LOS velocities on an absolute scale, our 
Doppler shifts have to be calibrated to high precision and cor- 
rected for the systematic effects introduced by the Earth's rota- 
tion, the orbital motion of Earth around the Sun, the solar rota- 
tion, as well as for the solar gravity redshift. 

The final accuracy and precision have to be high enough to 
be able to measure typical photospheric Doppler shifts which, 
in the case of the Si i line, correspond to velocities in the range 
of few hundred ms~'. Since our spectral range comprises two 
telluric JLO lines, we can use these to calculate our sampling 
(A per pixel). A Gaussian function with six terms was fit to 
the deepest part of these lines in our data in order to determine 
the position of their centers. The difference, Ax, is the distance 
between both lines in pixel units. The same process was car- 
ried out on the Four ier Transform Spectrometer (FTS) spectrum 
dKurucz et alJ [l984) from the Kitt Peak National Observatory, 
obtaining AA = 1 .873 A for the distance between both lines. The 
spectral sampling is merely calculated by dividing AA/Ax. This 
was done for each map using the average spectrum from a small 
non-magnetic region (containing about 1000 pixels). The mean 
spatial sampling of all maps is 1 1.035+0.010mA/px. Ideally, the 
sampling should always be the same. It can be theoretically cal- 
culated using the specifications of the telescope and of the spec- 
trograph. A comparison with the theoretical sampling (11.240 
mA/px) reveals that the calibration delivered a slightly smaller 
value, but not too different. 

We are now able to construct the wavelength array using 
the calculated spectral sampling and a telluric line whose cen- 
tral wavelength is known, as a reference. According to the solar 
spectrum atlas of ISwensson et al.1 d 19701) . the telluric line clos- 
est to the He i triplet has a wavelength of 10832.120 A, although 
the authors also provide another value: 10832. 150 A. However, 
when we calculate the line center position using the Gaussian 
fit to the FTS spectrum we find a wavelength of 10832.099 A. 
Note that the wavelengths of the FTS atlas are not corrected for 
the gravity shift. This is perfectly adequate for telluric lines, 
which are, indeed, not affected by it. Other wavelength val- 
ues, differing slightly from those mentioned previously, can also 
be found in the literature. These discrepancies led us to make 
our own estimate for the center wavelength of the H2O telluric 
line. The process, explained in Appendix [B] results in a wave- 
le ngth of 10832. 108 A This value is supported by the work 
of lBre ckinridg e"& Halll (Il973l) . who inferred a wavelength of 
10832.109 A with an accuracy approaching +1 mA. 

The newly constructed wavelength array is referred to a ter- 
restrial reference frame, that does not account for any relative 
orbital motions. We foll owed the calibration proced ures pre- 
sented in Appendix A of M artinez Pillet et al] (1 19971) . and ref- 
erences therein, adapted to the Observatorio del Teide, to obtain 
absolute line-of-sight velocities. Line shifts due to Earth's rota- 
tion, orbital motion of Earth around the Sun and solar rotation 
have been corrected. In the same way, the gravitational redshift, 
AAq = {GM q IRqC 2 )A (which translates into 23 mA for this spec- 
tral range), was also corrected in our calibration. The effect of 
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convective blueshift has been studied for the photospheric Si i 
line, however, after reviewing the literature and studying the re- 
sponse function to various physical perturbations of this line, we 
concluded that the correction is rather negligible owing to the 
formation height of the Si i line, which happens at a consider- 
able height (log r 2) above the surface. 

Appendix B: Determination of a new telluric line 
wavelength 

As mentioned in the previous appendix, the literature quotes 
several different values for the central wavelength (Aj) of the 
H2O telluric line next to the He 1 triplet. We calibrated our spec- 
trum using this telluric line as a reference, but the resulting ve- 
locity maps were clearly shifted to the red (when using Aj = 
10832.120 A) or to the blue (when using Aj = 10832.099 A). 
This is, we found systematic photospheric redshifts or blueshifts 
(depending on the wavelength used for the calibration) in the 
faculae, w here velocities are expecte d to be around zero (see 
Fig. 12 in iMartinez Pillet et alJ 1 199*71 where high filling factor 
faculae show no velocity shift). Shifts of the same order of mag- 
nitude were found when compared to new data from a recent ob- 
serving campaign in 2010 with the Tenerife Infrared Polarimeter 
at the VTT. We attribute this inconsistency to an incorrect value 
of the central wavelength of the telluric line and therefore we 
corrected its wavelength using the following method: we took 
three fiat fields from our 2010 campaign with TIP-II (August 
21 and 22); two from the morning and one from the afternoon. 
All of the fiat fields were taken at disk center in the quiet sun, 
with a random circular movement of the telescope pointing of 
up to 50". These maps were then used as input data for the stan- 
dard reduction procedure including flat field, dark curren t and 
polarimetric calibration corrections (Collados 1999, 2003). The 
spectral sampling for each flat field was inferred using the same 
procedure described in Appendix lAl Using Aj = 10832.120 A 
as the reference telluric line, a mean redshift of the Si 1 line of 
AA ~ 0.01 17 A was found. As mentioned above, this line is not 
expected to show any convective blueshift due to its formation 
height. Also no redshifts are expected at disk center. Thus, this 
systematic redshift was subtracted from the Aj = 10832.120 A 
telluric line, yielding a new wavelength of A^ EV/ = 10832.108 

A. With the new reference, the average photospheric facular ve- 
locity for all maps resulted in 0.06 kms -1 . Since the faculae 

observed here are very compact and have large filling factors 
(typically, higher than 50%), this value is expected to be zero. 
We thus conclude that the systematic effects of our velocity cal- 
ibration are smaller than 60 m s _I or around 2 mA. 
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